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Abstract 

The fraction r(t) of spins which have never flipped up to time t is studied 
within a linear diffusion approximation to phase ordering. Numerical simu- 
lations show that, even in this simple context, r(t) decays with time like a 
power-law with a non-trival exponent 9 which depends on the space dimen- 
sion. The local dynamics at a given point is a special case of a stationary 
gaussian process of known correlation function and the exponent 9 is shown 
to be determined by the asymptotic behavior of the probability distribution of 
intervals between consecutive zero-crossings of this process. An approximate 
way of computing this distribution is proposed, by taking the lengths of the 
intervals between successive zero-crossings as independent random variables. 
The approximation gives values of the exponent 9 in close agreement with the 

results of simulations. 
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Phase-ordering and domain growth in systems quenched from a disordered phase to 
a two-phase coexistence region has been for a long time a subject of active research [|TJ. 
Recently, a new facet of this problem has been uncovered 0-0- For an Ising model initially 
in a random configuration, and evolving according to zero temperature Glauber dynamics, 
the fraction r(t) of spins which have never flipped up to time t decreases like a power-law 
II 

r(t)~t- 9 . (1) 

Numerical data indicate that the exponent 9 is non-trivial and varies with the dimension of 
space @||. Similar exponents have been found in reaction-diffusion models with particles 
of two different mobilities ||[7||| . The decay ([!]) has also been observed for the fraction of 
surface that has never been wet in an experiment on the growth of breath figures |J. 

The simplicity of Glauber dynamics and its relation to soluble voter models has made 



possible an exact determination of 9 for Id Ising and Potts models JTO] . In higher dimensions 



only numerical simulations P,|TT] or approximate methods |T2j have been used so far. A very 
simple (and rather successful) description of coarsening for a non-conserved order parameter 
was proposed some time ago by Ohta, Jasnow and Kawasaki (OJK) ||13|| . In this OJK 
approximation, the order parameter configuration <p(f, t) a time t after the quench, is the sign 
of a gaussian field, initially random, which evolves according to a linear diffusion equation. 
Namely, 

<f>{f,t) = siga[A{f,t)] (2) 

with 

d t A(f, t) = V 2 A(r, t) and (A(f, Q)A(P, 0)) = 5 d (f - ?) (3) 

In this letter, we first report numerical simulations showing that even within this approx- 
imation, the fraction r(t) of persistent spins decays (|1|) with a non-trivial 9 which varies with 
the spatial dimension d. We then show that determining r(t) for the OJK dynamics is a par- 
ticular case of a first-passage problem for correlated random variables. In a logarithmic time 



scale, these problems can be mapped onto stationary gaussian processes. The exponent 9 is 
given by the asymptotic decay of the distribution of intervals between zero-crossings for the 
associated process. To our surprise, given the correlation function of the gaussian process, 



the determination of this asymptotic decay turns out to be a hard unsolved problem |T^-|T6 
We propose here a simple approximation which consists in taking the lengths of successive 
intervals between zero-crossing as independent identically distributed random variables. The 
probability distribution of the interval lengths is then chosen so as to reproduce the known 



sign-correlation function of the stationary gaussian process [[17]]. 

We have simulated the diffusion equation (|3|) by using a finite difference first-order for- 
ward Euler scheme on cubic lattices in d = 1, 2, 3 



Ai(t + At) = (1 - 2dAt)Ai(t) + AtY, M*) ■ ( 4 ) 

<ij> 

We used mostly At = 0.1 on lattices of up to 10 6 sites and we measured r(t), the fraction of 
lattice sites i for which Ai never changed sign up to time t. Our results are shown on Fig.l 
and indicate that r(t) decays as in ([]]) with 9 ~ .12 in d = 1, 9 ~ .18 in d = 2 and 9 ~ .23 
in d — 3. Other choices of At indicate that 9 is not affected by the value of At (as long as 
At < l/2d to ensure the stability of the Euler scheme). 

The stochastic process at a particular lattice point (the origin, say) is of the type 

A(t) = y £K(x,t)r }x (5) 

where the rj x are random gaussian variables (with (77^) = 1) and the kernel K(x,t) is 

K(x, t) = Const x^exp(-x 2 /4t)/(47rt) d/2 . (6) 

Other examples with different kernels include usual random walks (A = rjt) for which K is 
the Heavyside function K(x,t) = H(t — x) and 9 = 1/2, and walks under the influence of 
a random force (A = rjt) for which K(x,t) = (t — x)H(t — x) and 9 = 1/4 ||i8| , [T9"|| . Clearly, 
the first passage exponent 9 depends on the shape of the kernel K(x, t) and takes differents 
values for square, triangular or diffusive kernels. 



It is useful to note that (]5|) is a gaussian process since it is a linear combination of random 
gaussian variables. This process is therefore entirely characterized by its pair correlation 



function (A(t)A(t')). This is also true of the rescaled variables X(t) = A(t)/^(A(t)A(t)) 
which can be considered as well since we are only interested in the sign changes of A(t). 
In the long time limit, the correlation function (X(t)X(t')) depends only on the time ratio 
t/t' for the process considered (J|). In a logarithmic time scale, u = \og(t), if one defines 
Y{u) = X(t), Y{u) is therefore a gaussian stationary process 

{Y(u)Y(u')) = C(u - u) (7) 

with C(u) = [1/ cosh(w/2)] d / 2 for the OJK process in dimension d (and C(u) = exp(— \u\/2) 
for the usual random walk). C is a rapidly decreasing function of u. Therefore, Y(u) and 
Y(u + l) are effectively uncorrelated for sufficiently large I and it is intuitively clear that the 
probability that Y(u) keeps the same sign for an interval of length I decreases exponentially 
with / for large I 

Prob{Y(u + I') > 0, < I' < 1} ~ exp(-ei) for / > 1 (8) 

and 9 measures the (inverse) decorrelation time in the w-variable. The asymptotic behavior 
(§) implies that the probability of not flipping in the original t variable from time t\ to 
with t\ <C ti is a power law (t\jt2) 6 as numerically observed above. 

In order to compute 6, one needs the large I behavior of the distribution p(l) of the lengths 
I of intervals between consecutive zeros of a gaussian stationary process characterized by a 
given C{1). As observed by Rice more than forty years ago [fL4}| , this is a surprisingly hard 
problem which is still unsolved |T5|JT8|JT6|1 . So, one has to resort to approximate methods. 
Perturbative and variational techniques using the free brownian walk as a starting point have 
been proposed recently [12|,20l. However, the free brownian walk and other similar markovian 



processes have an infinite density of zero-crossings and are not a good starting point for cases, 
for which C"(0) is finite like the OJK process, which have a finite density of zero-crossings 



(the density of zeros is given [14] by J—C"(0)/n). We propose here an approximation more 



suited to these cases which is inspired by recent works on the distribution of domain sizes in 
a Id coarsening Ising model |^TJ where a similar approximation has been shown to be rather 



accurate ||22|| . It consists in considering the lengths between consecutive zero-crossings as 
independent random variables, identically distributed according to a probability distribution 
p(l). This p(l) is determined by requiring that the sign correlation function is the same for 
the approximate independent interval approximation and for the original stationary gaussian 
process p3 . 



For a process with correlation function C(l), the probability P+(l) that Y(u)Y(u + l) > 
is given by []14 | 



P+(l) = - + -arcsin(C(/)) . (9) 
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To obtain the corresponding probability for the independent interval process, it is convenient 
to determine first the probability P cr (l) that a point at distance / from an up-crossing be 
positive. Either the first interval after the up-crossing is longer than I, or the interval is of 
size l\ < I and the point is positive with probability (1 — P cr (l — h))- So P cr (/) satisfies 

/+oo A 
dh p(h) + / dh p(h) (1 - P cr (l - h)) . (10) 
Jo 

Taking Laplace transforms (which we denote with tilde), one obtains 

P cr (A) = [A(l+p(A))]- 1 . (11) 

This can now be used to determine the probability that Y (u)Y (u+l) > for the independent 
interval process. A given point is at a distance h from the first crossing on its right, with 
probability 

nn , ti°°p(l) dl - p(A) - 1 . 

If Y(u) and Y(u + I) have the same sign either they belong to the same interval or Y(u + I) 
has a sign opposite to points immediately to the right of the first crossing after u, so that 

/+oo A 
dh Q(h) + J o dh Q(h) [1 - P cr (l - h)} ■ (13) 



This gives for the Laplace transforms 

P+(\) = l/\-Q(\)P a .{\) . (14) 

Substituting flTT| ) and flT2"| ) into (P^), one finally relates the distribution p{l) of interval 

lengths to the probability P+(l) that two points at distance I have the same sign for the 

independent interval process (or rather their Laplace transforms) 

_ 1 - \\!P + (X) - A] g(0) 
P(A) - 1 + |VP + (A) - A] m ( ' 

For an integrable distribution of interval lengths, p(X) should tend to zero as A —>■ oo. Using 



(Hp, this gives p'(0) = 1/Pi.(0) which simply means that the mean interval size (— p'(0)) is 
the inverse of the density of zeros (— P( (0)). 

Equation ([TJ) would be exact if the successive intervals were uncorrelated. Here, the 
heart of the approximation is to choose the exact P+(l) given by (|9|) for a given process 
and determine p(l) through (pp . The large I behavior of p(l) is given by the singularity 
of its Laplace transform with the largest (negative) real part i.e. by the largest zero of 
P|(0) = A-A 2 P + (A) or, using (|), 

i r+oo 

y/-C»(p) = -9 / dlexp{9l) C\l) [1 - C 2 (l)}- 1/2 . (16) 

J 

In the OJK case, C(l) = [l/cosh(//2)] d / 2 . A numerical solution of @ gives 6 ~ .1203 for 
d — 1, ~ .186 for d = 2 and 6* ~ .2358 for d = 3 in very good agreement with the previous 
numerical results. For large d the independent interval approximation (|16]) gives 9 ~ k 
with k ~ .1455 ||24|| . Note that for the case of a free brownian motion, C"(0) = oo and the 
only possibility to satisfy (|T6|) is that the integral on the rhs diverges, leading to 9 = 1/2 
which is exact. 

The whole distribution p{l) can also be obtained by inverting the Laplace transform. 
One can obtain the Taylor expansion of p{l) around I = easily from (|I~5| ) and then draw 
p{l) using Pade approximants. The obtained shapes for the diffusive case in d— 1,2,3 are 
plotted in Fig. 2 and compared with results of direct simulations of the stationary gaussian 
processes. 



We have shown that non-trivial exponents for the number of persistent spins appear even 
in the simple OJK approximation of coarsening. Contrarily to the case of Glauber dynamics, 
9 increases with dimension. This means the OJK theory is too crude an approximation to 
predict r(t) for Glauber dynamics. It is remarkable that the independent interval approx- 
imation gives here such accurate predictions f° r the OJK exponents. Still, it would 
be interesting to better understand it, for example by a maximum entropy argument. Of 
course, it would be also interesting to estimate the correlations p(h,h) between successive 
intervals and see whether one could develop a systematic method of including these correla- 
tions. Another issue would be to see under what conditions, the distributions p{l) obtained 
via ( |15D are positive for general kernels K (as nothing guarantees a priori that it is so). 

Lastly, one could try to extend the OJK approximation to cases where the order pa- 
rameter is more complicated than a scalar (as in the Potts model) to see how the above 
aproximation can be adapted to such cases. 

Acknowledgments We are grateful to D. Dhar for many enlightening discussions on 
the subject of this paper when two of us (BD and VH) were visiting the Newton Institute 
(Cambridge, UK). We also thank A. Bray for informing us of his parallel work and J. P. 
Bouchaud for interesting comments. 



7 



REFERENCES 

[1] J. Langer in Solids far from Equilibrium, C. Godreche ed., (Cambridge U. Press, Cam- 
bridge, 1992); A.J. Bray, Adv. Phys. 43, 357 (1994) 

[2] B. Derrida, A. J. Bray and C. Godreche, J. Phys. A 27, L357 (1994) 

[3] D. Stauffer J. Phys. A 27, 5029 (1994) 

[4] A. J. Bray, B. Derrida and C. Godreche, Europhys. Lett. 27, 175 (1994) 

[5] P.L. Krapivsky, E. Ben-Naim and S. Redner, Phys. Rev. E 50, 2474 (1994) 

[6] P.L. Krapivsky, S. Redner and F. Leyvraz, Phys. Rev. E 51, 3977 (1995) 

[7] E. Ben-Naim, L. Frachebourg and P.L. Krapivsky Phys. Rev. E 53, 3078 (1996) and 
preprint 96 

[8] J. L. Cardy J. Phys. A 28, L19 (1995) 

[9] M. Marcosmartin, D. Beysens, J. P. Bouchaud, C. Godreche and I. Yekutieli, Physica 
A 214, 396 (1995) 

[10] B. Derrida, V. Hakim and V. Pasquier, Phys. Rev. Lett. 75, 751 (1995); J. Stat. Phys. 
(to appear) 

[11] B. Derrida, P.M.C. de Oliveira and D. Stauffer, Physica A 224, 604 (1996) 
[12] S. N. Majumdar and C. Sire, preprint (April 96) 

[13] T. Ohta, D. Jasnow and K. Kawasaki, Phys. Rev. Lett. 49, 1223 (1982) 

[14] S. O. Rice, Bell Syst. Tech. Journal 23 and 24 reprinted in N. Wax Selected Papers on 
Noise and Stochastic Processes (Dover, N. Y., 1954). 

[15] M. Kac, SIAM Rev. 4,1 (1962) 

[16] H. Cramer and M. R. Leadbetter Stationary and Related Stochastic Processes (Wiley, 

8 



1967) 

[17] After an oral presentation of our results at the 96 Nordic Stat. Phys. Workhop in Copen- 
hagen, A. Bray informed us that he and his collaborators have obtained independently 
results similar to ours: S. N. Majumdar, C. Sire, A. J. Bray and S. J. Cornell preprint 

(May 96) 

[18] H. P. McKean Jr., J. Math. Kyoto Univ. 2 227 (1963); M. Goldman, Ann. Math. Stat. 
42 2151 (1971) 

[19] Y. Sinai Theor. Math. Phys 90, 323 (1992), T.W. Burkhardt, J. Phys. A 26, L1157 
(1993) 

[20] V. Hakim and R. Zeitak (in preparation) 

[21] P. A. Alemany and D. ben-Avraham, Phys. Lett. A206, 18 (1995) 
[22] B. Derrida and R. Zeitak, preprint (April 96) 

[23] A. Bray has informed us that this approximation had been previously proposed for 
stationary gaussian processes by J. A. McFadden, IRE Transaction on Information 
Theory,IT-4, 14 (1957) 

[24] The constant k is solution of the equation 2\/2 = k J +oc duexp(K,u)u/ Jexp(u 2 /8) — I. 



9 



FIGURES 

FIG. 1. Log-log plot of the fraction r(t) of sites i for which has never changed sign vs. time 
for the discretized diffusion equation in dimension d = 1,2,3. r(t) shows a power-law decay with 
an exponent which depends on d. 

FIG. 2. The probability distribution function of having an interval of length I for the gaussian 
stationary process associated to the diffusion equation in d = 1, 2, 3 as obtained in the independent 
interval approximation (lines) and by direct simulations of the gaussian processes (d = 1 circles, 
d = 2 plusses, d = 3 diamonds) in Fourier space (1000 realizations of length 3276.8; bin size 0.1). 
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